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1 Introduction 

Quantum electrodynamics can be considered to be one of the most successful theoretical 
concepts, allowing to predict experimental results with unchallenged precision. Here the 
anomalous magnetic moment of the electron (see, [1] and references therein), positron- 
ium spectroscopy [2] or properties of light hydrogenic bound states [3] may be listed as 
most outstanding examples. At the same time the detailed investigation of its mathemat- 
ical structure has lead to fundamental field-theoretical concepts, like renormalization or 
the renormalization group with anomalous dimensions and the /3-function as important 
elements for the analysis of its high energy behaviour [4-6]. 

As a consequence of Ward identities the /3-function of QED is closely related to the 
vacuum polarization function H(q 2 ) and, in contrast to non-abelian gauge theories like 
QCD, no vertex corrections need to be evaluated. Indeed, the coefficient of the n-th term 
of the perturbative series of the /3-function can be directly obtained from the absorptive 
part of n(q 2 ) at the same order. The same absorptive part, evaluated including QCD 
corrections, can be employed to predict the Adler function [7] and thus the total cross 
section for electron positron annihilation into hadrons or the r-lepton semileptonic decay 
rate. The precise experimental information for these fundamental quantities has stimulated 
dedicated projects for the evaluation of three- [8], four- [9, 10] and recently even five- loop 
corrections [11-15] for these fundamental quantities. Conversely, it is obvious that the 
information contained in these five-loop QCD results is sufficient to evaluate the QED 
vacuum polarization and the f3 function in the same five-loop order. Partial results in this 
direction have been presented in [11, 13, 16]. 
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It is the purpose of this work to collect all available results for the vector correlator in 
massless QCD in a form suitable for a complete evaluation of the QED f3 function. Two 
cases will be considered: the QCD contribution to the QED /3-function of order a 2 a^ and 
the QED ^-function of order a 6 for a generic number of charged fermions. 

For readers's convenience we provide in the paper the complete results for the vector 
current correlator, related /3-functions and anomalous dimensions, including also lower 
order contributions. As the latter are known since long, we refer the reader to works [17— 
19] for corresponding historical discussions and citations. Our full results are also available 
(in computer-readable form) in http://www-ttp.physik.uni-karlsruhe.de/Progdata 
/ttpl2/ttpl2-018. 

2 Preliminaries 

Our main object will be the polarization function H(L, a s ) of the flavor singlet vector 
current defined as 

{~g a pq 2 + q a qp) U(L, a s ) = i f d^xe 1 ** {0\Tj a (x)jp(0)\0), (2.1) 



with j a = i>ilai>i and Q 2 = -q" 2 



<*s(M) 



We work within massless QCD 



with rif quark flavours, \i is the normalization point of the MS-scheme. 




Figure 1: Examples of a lowest order and five- loop non-singlet (a),(b) and singlet (c),(d) 
diagrams contributing to the vector correlator. 



All diagrams contributing to H(L,a s ) can be naturally decomposed into two classes: 
singlet and non-singlet ones. The non-singlet diagrams are those with both external cur- 
rents belonging to one and the same quark loop; all other diagrams are referred to as singlet 
ones (see Fig. 1). As a consequence of this classification the full polarization function can 
be presented as follows: 



U(L,a s 



(2.2) 



n f U NS (L,a s ) + nj U SI (L,a s ). 

It is worth to note that the knowledge of II^ 5 and II s7 is enough to construct the polar- 
ization function arising from two vector currents with generic flavour structure. Indeed, 
if j% = 4>lat a ip an d ja = ip"f a t b ^ (i° and t b being two arbitrary matrixes in flavor space) 
then the corresponding polarization function U ab is, obviously, given by the relation 
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For instance, the polarization function corresponding to the electromagnetic vector current 



■EM 

J u 



^2 QtiPa^i 



(qi stands for the electric charge of the quark field ipi) is given by the expression: 

2 

(2.3) 



The absorptive part of H EM is directly related to the important physical quantity R(s) = 
a(e + e~ — > hadrons)/cr(e + e _ — > via the relation 

R(s) = 6vr (U EM (-s - ie) - U EM (-s + ie)) . (2.4) 

Another useful object - so-called Adler function — is defined as [7] 

d n 2 

D(L,a s ) = -12vr 2 Q 2 — ^U(L = ln^,a fl ). (2.5) 



Note that the choice of the coefficients 67r and — 127r 2 in eqs. (2.4) and (2.5) is conventionally 
fixed by the requirement that in Born approximation 

R(s) = (^2 ?n d R C 1 + °^)) . D ( L ' °«) = n f d R i 1 + °(«-)) > 

where the d^j is the dimension of the quark representation of the colour gauge group (djt = 3 
for QCD). 

In spite of the fact, that the vector current in QCD is a scale-invariant object (that is it 
has zero anomalous dimension), the polarization function is not due to the short distance 
singularities of the T-product in eq. (2.1). The corresponding evolution equation reads 
(see, e.g. [17]) 



2 d . „, s d 



^ 9^ + /3(as Wj n(L ' as) = 7(as) (2 ' 6) 
or, equivalently, 

d d 

— U(L, a s ) = 7(a fl ) - /9(a s )— IL(L, a s ). (2.7) 
oL oa s 

Here the anomalous dimension 7 is a series in a s of the form 



16 vr 2 

i>0 



and the QCD /3-function is defined as 

,_d 



^A_a s = P{a s ) = -Y J ^ s +2 - (2.9) 



For our analysis we need to know the QCD /3-function with three-loop accuracy, the corre- 
sponding result is known since long from [20, 21]. Note that the ratio R{s) and the Adler 
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function are scale-invariant as the rhs of eq. (2.6), being an anomalous dimension does not 
depend on either Q, \i or L. 

In massless QCD eq. (2.7) directly leads to an important relation for the Adler function: 

D(L, a s ) = 12vr 2 U(a s ) - (^3(a s )-^j U(L, a s )^j . (2.10) 

Eqs. (2.7) and (2.10), are crucial for our ability to compute the Adler function and R(s) 
at order oQ which, in principle, are determined by five-loop diagrams. The simplification 
arises because the rhs of eq. (2.10) is, in fact, expressible exclusively through four-loop mass- 
less propagators. First, only the four-loop O(a^) approximation to IT(L,a s ) is required as 
the /3-function starts from order a 2 . Second, the evaluation of any (L + l)-loop anoma- 
lous dimension in the MS-scheme can be reduced, with the help of the i?*-operation [22], 
to the evaluation of some L-loop masslesss propagators [23]. Finally, four- loop massless 
propagators can be reduced to 28 master integrals. The reduction is based on evaluating 
sufficiently many terms of the 1/D expansion [24] of the corresponding coefficient functions 
[25]. The master integrals are known analytically from [26, 27]. 

Needless to say that any direct way of computing 5-loop diagrams, contributing, for 
instance, to the Adler function without the use of (2.10) is hopelessly difficult and, presum- 
ably, will stay so for years and years ahead. 



3 Results for the vector correlator 

In this section we present all currently available results for the polarization function II 
and the anomalous dimension 7 which have been utilized in works [12-15] to construct the 
Adler function and the ratio R(s) in massless QCD to order a^. 

3.1 The polarization function II 

By presenting the perturbative expansion of the polarization function II as follows 



IT 



dR 



16vr 2 ^ 

i>0 



IT 



16vr 2 




(3.1) 



we get 



r> NS 

Po 



P? S 



n NS 
P2 



20 

y 

C F 



Cl 



12 " 4C3 



143 37 A ' 
-W— 6 6 + 106 



+Cf Ca 



44215 
2592 



— C - V 
18 ^ 3 W 



+ C F Tn 



3701 38 , 

\ Cs 

648 9 s 



(3.2) 
(3.3) 

(3.4) 
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P3 



c% 



31 13 245 . ' 

-i92 + y C3 + — Cs " 35C7 



+T 2 n 2 f C F 



+Tn f C 2 F 
+T rif Cp C a 
+C F Ca 



7505 1553 > 



11 



196513 
23328 
250 



809 _ 20 ' 
loir 3 ~ ~9 



10368 



54 



4(3 + ^4-^5 



5559937 41575 
+ 



93312 
382033 46219 



/- j /-^ 

1296 W 3 ^ 3 



11 515 
24 C4 + "27~ C5 



+Cf C a 



r> SI 



20736 
34499767 _ 
373248 

d abc d abc t 431 



Ca 



864 
147473 



11 9305 35 
^C4 + W C5 + yC7 



55 2 11 
2592~ C3 + "6" C3 + 48 C4 



28295 
864 



-Cs 



35 
12 



21 



1 



1 



64 Cs "6 C3 " 16 C4+ 16 C5 



(3.5) 



(3.6) 



d R { 1728 

Here Cp and Ca are the quadratic Casimir operators of the fundamental and the adjoint 
representation of the colour Lie algebra, d abc = 2Tr({^-, ^}4r), T is the trace normal- 
ization of the fundamental representation. The exact definitions of the colour structures 
^abcdjabcd an j jabcd^abcd a p pear i n g below are given in [28]. For QCD (colour gauge group 
SU(3)): 



d R = 3 , C F = 4/3 , C A = 3 , T = 1/2 , 



^abcd ^abcd 



15 5 40 

iabcd jabcd " jabc mbc u 

Y> d F d F -^2' d d -y 



(3.7) 



For the particular case of the U(l) gauge group the colour factors assume the following 
values: 



d R . = 1 , C F = 1 , C A = , T = 1 , df cd dT = , df cd df cd = 1, d abc d abc = 16. (3.8) 

Note that in eqs. (3.2-3.6) we have set to zero L = ln(/^ 2 /<5 2 ). The full dependence 
on L can be easily restored from evolution eq. (2.7) and the anomalous dimension 7 given 
below. 

3.2 The anomalous dimension 7 

On decomposing the anomalous dimension 7 into non-singlet and singlet terms 



7 = n /7 " s + nJ 7 s ', 7^ S 



dR 
16vr 2 



i>0 



16vr 2 



i>3 



i a s I 



(3.9) 



we get 



~,ns 
7o 



3' 

71 = 

72 — °F 



133" 
144 



+C F Tn 



I 



11 

36 



(3.10) 
(3.11) 
(3.12) 
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^NS _ f~<3 

73 — U F 


23" 
~32 


+T 2 n}C F 


77 ' 
972. 


+Tn 




169 11 A 

1 C3 

216 9 




+Tn f C F C A 




769 


-g-Cs 


+c| 


Ca 


"215 11 






3888 


216 ~ 18 



+C F C 2 a 



5815 11. 

1 C3 

15552 18 s 



(3.13) 



C% 



4157 1 

h -C3 

1536 2 Sd 



+n f 



jjabcd jabcd 
F F 



da 



13 4 10 

"12 ~ 3^ 3 + l"^ 5 



4 3 W 3 s5 



+T 3 n 3 f C F 



107 1 
7776 + 54 s3 



+ Tn f C% 



4961 
10368 
2509 



119 11 

108 + 24 ^ 



+T 2 n) C F C A 



8191 



1152 



67 145 
24 ~ ~24~ 



31104 
66451 2263 



563 _ 11 
432^ 3 ~ 24 ^ 



+ Tn f C F C 2 A 



22423 9425 143 15 
31104 ~ 1728^ + ~9() 4 + 32^ 5 



13824 
+C f Ca 



364 s3 96 C4+ 16 < ' 5 



2585 
256 



71 935 ' 
12 s3 + ~96~^ 5 



+ C 2 F C 2 A 



+ c F c A 



882893 11501 121 715 
82944 + 3456 ^ 3 + 192^ 4 ~ ~6A < " 5 



1192475 5609 
497664 + 3456^ 



121 275 
192^ + 128^ 5 



(3.14) 



73" 



7f 



d abc d abc J u 

1 144 



6 



(3.15) 



C F 



13 _ 1 5 " 
"48 ~ 3^ 3 + f) 



C A 



1015 
2304 



659 , 11 5/ 
768 C3 + M C4 + 64 C5 



+Tn f 



55 41 

1 C3 

576 192 s 



1 

— ( 
16 



48 



Cs 



(3.16) 



4 QED /3-functions in five-loop order for different schemes 



4.1 Massless QCD and QED: MS-scheme 

The polarization function 14 is known to be directly related to the photon propagator, 
namely 1 

9 d(Q 2 ) 1 

Vafi(q) = {g^q 2 -q a q/3) d(Q 2 ) = - - ^ ngM(Q2) - (4-1) 

The combination e 2 d(Q 2 ) is often referred to as "invariant" charge as it is renormalization 
scale and scheme independent due to the corresponding Ward identity. The independence 



1 Without loss of generality we will use the Landau gauge for the photon field. 
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of the invariant charge on the renormalization scale fi directly leads to the RG equation 
for the QED coupling constant A(fi) = a(/i)/(47r) = e(/x) 2 /(167r 2 ): 



2 d 



A = f3 EM (A, a s ) = A 2 (16vr 2 ) 7 BM (a s ) = A 2 d R ^ 7, 



(4.2) 



i>0 



with 



7 EM (a s ) 



E ft 2 ) ^ NS + (e *) 2 ^ Tf M - (E ft 2 ) tf s + (e ft) 2 ^ 



The /3-function /3 BM (A, a,,) describes the QCD-induced corrections to the running of a in 
the MS-scheme. 

Using now eqs. (3.10-3.16) and substituting the values of the colour factors correspond- 
ing to the SU(3) colour group we find: 



^ 1 (A,a s ) = A 2 (^ 2 ) |4 + 4a s + a 2 



o / 10487 110 
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216 ~ ~27 



<3 
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~ 18 
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~ 972 
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4 / 2665349 182335 605 ^ 31375 
F ° s \ 41472 + ~86^ C3 " T6" C4 " ^88~ C5 

11785 58625 



+ n s 



+ n 2 



715 13325 
864 ^ 3 + l^P 4 + 432 ^ 
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4729 3163 55 
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+ n 



f 
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)*l 









55 20 



54 9 ^ J / ' s \ 864 



+ n f 



<3 +< 

275 205 



15552 108 
11065 34775 



864 
25 



<3 



+ 



55 3875 
¥ C4 + W Cs 



f 144 ^ 3 12 ^ 36^ 



(4.3) 



For the particular cases of 4, 5 and 6 quark flavours eq. (4.3) takes the form (the normal- 
ization is chosen in such a way to facilitate the comparison with [19]) 

/3 BM (A,a s )\ nf=4: = A A 2 (l.lll + 1.111a, + 2.214a 2 + 1.210a 3 - 3.904a^) , (4.4) 

f3 EM (A,a s )\ nf=5 = 4 A 2 (1.222 + 1.222 a s + 2.249 a 2 - 1.227a 3 - 13.429a 4 ) , (4.5) 

f3 EM (A,a s )\ nf=6 = 4 A 2 (1.667 + 1.667 a s + 2.813a 2 - 5.791a 3 - 32.336af) . (4.6) 

Another case of interest is pure QED, that is a theory with rif single-charged fermions 
minimally coupled to the photon field. In this case the corresponding EM current is 
identical to the flavour singlet current. The corresponding f3- function is, obviously, obtained 
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from the general formula (4.2) by taking the QED values for the colour factors and setting 
qi = 1 and a s (fi) = 4i(/j): 



^ ED {A) = n f 



A A 2 



+ A 5 
+ A 6 



n 



+ n 



760 2 
-4Sn f + — n f - 

+ 128C 3 
16000 



+4rifA 
832 



A 4 



44 



2n / + lT n / 



9 



Can} 



1232 
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n 3 

n f 



4157 



6 

21758 



+ n1 



7462 



992C 3 + 2720C 5 



81 



. 416^ 



1280 



856 128 , 

1 

243 27 s 



(4.7) 



If we set nt = 1, then the above result takes the form: 



P QED {A) 



-A 2 + 4A d 



62 



5570 832 



243 



'<2 



+ A b 



195067 800 416 6880 ' 
486 JT'' 3 ~ + 3 ^ 5 



(4.8) 



or, numerically, 

= -A 2 + 4 A 3 
3 

= -A 2 + 4 A 3 



62 

yi 4 +i 5 (-7.116 - 126.93 s/ ) + A 6 (776.39 + 729.63 s/ ) 
62 

—A 4 - 134.045 A 5 + 1506.02 A 6 , 



(4.9) 



where in the first line we have explicitly separated non-singlet from singlet contributions. 

For future reference it is also useful to present the evolution equation for the QED 
polarization function: 



— U QED (L, A) = 7 QED (A) - f3 QED (A)-^-n QED (L, A). 
oL oA 



Here 



f3 QED (A) = A 2 (16tt 2 )j qed (A), 

w* D (L, A) + n 2 f U SI ' QED (L, A). 



U QED (L,A) = n f U 



(4.10) 



(4.11) 
(4.12) 



In addition, i[ NS ' QED , U SI ' QED and ^ QED are IP^, Il s/ and 7 respectively with all colour 
factors substituted according to eq. (3.8) and a s = A A. 

4.2 Massless QCD and QED: MOM-scheme 

The MOM-scheme for the QED coupling constant is defined by the requirement that at 
Q 2 = [i 2 the polarization function would vanish. The scheme independence of the invariant 
charge directly leads to the following relation between MS-renormalized QED coupling 
constant A(fj.) and its MOM analog 



A(jm) 



1 + {Att) 2 A(fi)U(L = 0,a a (fi))' 



(4.13) 
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The above equation allows to express easily the QED /3-function in the MOM-scheme 
via the corresponding MS /3-function and the perturbative expansion of the polarization 
function as given in eqs. (2.3,3.2-3.6,4.3). After differentiating the rhs of (4.13) with respect 
to fx and using eq. (2.6) we obtain 



(3 EM (A,a s ) = 16tt 2 A 2 
4 



7 EM (a s )-P(a s ) ^-U EM (L = 0,a s 
da. 



where 



-A 2 D EM (L = 0,a s ), 



D EM (L, a.) = -12 vr 2 Q 2 ^^(L = In ^, a,). 



(4.14) 
(4.15) 



Explicitly, for the SU (3) colour group we get 
A A 2 



{J EM (A,a s ) 
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(4.16) 



For particular cases of 4, 5 and 6 quark flavours eq. (4.16) reads 

f3 EM (A, a s )| n/=4 = 4 A 2 (1.111 + 1.111 a s + 1.694a 2 + 2.881 a 3 + 28.132a*) , (4.17) 

f3 EM (A,a s )\ nf=5 = 4 1 2 (1.222 + 1.222 a, + 1.723a 2 - 0.879a 3 + 10.703 a*) , (4.18) 

f3 EM (A, a s )\ nf=6 = 4 1 2 (1.667 + 1.667a, + 2.156 a 2 - 6.995 a 3 - 13.994 af) . (4.19) 

In the same way one could derive the pure QED /3-function (3 QED (A) in the MOM- 
scheme (it was first introduced by Gell-Mann and Low in [5] under the name u ift function" ) . 
Indeed, an analog of eq. (4.13) now assumes the following form: 

A<ji) 



1 + (4tt) 2 A(ji) Ii QED {L = 0, A) 



(4.20) 
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The resulting formula for the QED /3-function in the MOM-scheme reads (we have used 
relation (4.10) and definitions (4.11,4.12)) 



pQ ED (A) = 16vr 2 i 2 
4 



-A 2 D QED {L = 0,A), 
3 



(4.21) 



where D QED (L, A) = -12vr 2 Q 2 gg, I1 QED {L = In £, A). In addition, the coupling constant 
A(fi) appearing in the rhs of eq. (4.21) should be expressed in terms of A(fj,) using an 
inversion of eq. (4.20). Explicitly, for the U(l) gauge group we arrive at the following 
result: 
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(4.22) 



+ 4 



27064 26240 
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-C3 + 1024C3 2 + 2560C 5 



8576 3584 , 5120 , ' 



After setting nj = 1 we get: 



j8 5BD (A) = -i 2 + 4i 3 + i 4 
3 
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10 „ 256 . 1280 . 

186 + _ {3 __ C5 



122387 10112, t a2 23680 A 
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-,4 2 + 4 A 3 + 3.199 ,4 4 +A 5 (-26.918 - 126.929 s1 ) + A 6 (1054.41 + 413.592 s7 ) 
3 



A 6 
4 



-A 2 + 4 A 6 + 3.199 ,4 4 - 153.847.4 5 + 1467.998 A b 



(4.23) 



where the singlet contribution has again been identified separately. 
5 Discussion 



£ 


1 


2 


CO 


4 






Cs 


C3> C4> C5 


C3 , C4 , C5 , Cf ' C6 , C7 



Table 1: Possible irrational structures which are allowed to appear in £-loop massless 
propagators. 
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Let us discuss the structure of trancendentalities appearing in our results. It follows 
from work [26] that the variety of ^-constants entering into the MS-renormalized (euclidian) 
massless propagators depends on the loop order according to Table 1. Table 2 provides 
the same information about possible irrational numbers which could show up in anomalous 
dimensions. An examination of eqs. (3.2-3.6, 3.10-3.16) immediately reveals that the 



t 


1,2 


3 


4 


5 






Cs 


C3> C4, C5 


C3, C4, C5, Cf) C6) C7 



Table 2: Possible irrational structures which are allowed to appear in £-loop anomalous 
dimensions and /3-functions. 

observed pattern of trancendentalities is significantly more limited than what is allowed by 
Tables 1 and 2. Indeed, the four- loop anomalous dimension 73 contains no £4 and no £5 while 
the three- loop polarization function contains £5 but does not comprise £4. Let us move up 
one loop. The situation is getting even more puzzling: the five- loop anomalous dimension 74 
does contain £4 but still does not include £|, Cq and Qj. The four- loop polarization function 
contains £4 but is free from ^. Unfortunately, we are not aware about the existence of any 
reason behind these remarkable facts, except for one observation, namely the absence of £4 
in the MOM /3-functions (4.16) and (4.22). 

Indeed, according to eqs. (4.14) and (4.21) the constant £4 does not appear in these 
two /3-functions since it does not appear in the Adler function. However, the puzzle of the 
absence of £4 in 0{pq^) contribution to the Adler function has been recently fully explained 
in [26]. The explanation is based on a quite peculiar structure of irrational contributions 
to each four-loop master integral. 

Why this absence continues to hold at five loops is still a mystery (at least for us). In 
all probability it is connected with some regularities of five-loop master integrals. But here 
starts terra incognita . . . 

6 Conclusions 

We have presented four new results, namely the QED /3-functions /3 BM and f3 QED in the 
MS- and MOM-schemes. 

We have described the status of results for the vector correlator in massless QCD. 
These are not completely new as they have been used to produce the Adler function and 
R(s) in works [12-15]. Nevertheless, we believe that the separate presentation of the 
polarization function and its anomalous dimension is both useful and instructive. First, it 
reflects the real way how the calculations have been done. Second, it clearly demonstrates 
puzzling regularities of the structure of irrational terms contributing to II and 7. Third, 
it makes trivial the construction of the QED /3-function in the MS- and MOM-schemes 
(cmp. with the somewhat unnecessary complicated "inverse engineering" employed in [19] 
to reconstruct the function f3 EM (A, a s ) at four loop level). 
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The calculations of II and 7 have been performed on a SGI ALTIX 24-node IB- 
interconnected cluster of 8-cores Xeon computers using parallel MPI-based [29] as well 
as thread-based [30] versions of FORM [31]. For the evaluation of colour factors we have 
used the FORM program COLOR [32]. The diagrams have been generated with QGRAF 
[33]. The figures have been drawn with the help of Axodraw [34] and JaxoDraw [35]. 

This work was supported by the Deutsche Forschungsgemeinschaft in the Sonder- 
forschungsbereich/Transregio SFB/TR-9 "Computational Particle Physics" and by RFBR 
grants 11-02-01196, 10-02-00525. 

Finally, we want to note that the result (4.8) for the function /3 QED (A) with rif = 1 
was first reported by one of the present authors in September 2011 during the 10-th Inter- 
national Symposium on Radiative Corrections, RADCOR 2011 (see the 10-th page of the 
file http : //www. icts . res . in/media/uploads/Program/Files/chet . pdf ). 
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